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Abstract 

We use linearized time-dependent (current) density functional theory to study 
collective transitions between the two lowest subbands in GaAs/AlGaAs quan- 
tum wells. We focus on two particular systems, for both of which experimen- 
tal results are available: a wide single square well, and a narrow asymmetric 
double quantum well. The aim is to calculate the frequency and linewidth of 
collective electronic modes damped through electron-electron interaction only. 
Since Landau damping, i.e. creation of single electron-hole pairs, is not ef- 
fective here, the dominant damping mechanism involves dynamical exchange- 
correlation effects such as multipair production. To capture these effects, 
one has to go beyond the widely used adiabatic local density approximation 
(ALDA) and include retardation effects. We perform a comparative study of 
two approaches which fall in this category. The first one is the dynamical 
extension of the ALDA by Gross and Kohn. The second one is a more recent 
approach which treats exchange and correlation beyond the ALDA as vis- 
coelastic stresses in the electron liquid. It is found that the former method is 
more robust: it performs similarly for strongly different degrees of collectivity 
of the electronic motion. Results for single and double quantum wells com- 
pare reasonably to experiment, with a tendency towards overdamping. The 
viscoelastic approach, on the other hand, is superior for systems where the 
electron dynamics is predominantly collective, but breaks down if the local 
velocity field is too rapidly varying, which is the case for single-electron-like 
behavior such as tunneling through a potential barrier. 
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I. INTRODUCTION 



In recent years, time-dependent density functional theory (TDDFT)tJo has become the 
method of choice for describing dynamical electronic properties of atoms, molecules and 
solids in the linear regime and beyond. Until now, TDDFT has been mostly applied within 
the adiabatic local density approximation (ALDA)i for the dynamical exchange-correlation 
(xc) potential. However, some propositions have been put forward for improving upon the 
ALDA. The objective of earlier attempts@'i was to obtain approximations for the xc potential 
which would still be local in space, but not in time. All these approximations were found to 
suffer from inconsistencies, such as the failure to satisfy the so called "harmonic potential 
theorem" (HPT)H, and other basic symmetries. Only recently it has become clear that the 
root of these difficulties lies in the fact that the xc potential in TDDFT is an intrinsically 
nonlocal functional of the density, that is, a functional that does not admit an expansion 
in terms of the density gradients!'!. Consequently, the most recent extensions of TDDFT 
beyond the ALDA are formulated in terms of the local current density0 or the motion of 
local fluid elements of the electron liquidS. 

Among the variety of phenomena which have been investigated using TDDFT methods 
(see Ref. |]for a recent review), collective electronic excitations in semiconductor heterostruc- 
tures are currently of great theoretical and experimental interest. In particular, much effort 
has gone into the study of parabolic quantum wells with and without imperfections@~0. 
However, the only collective mode that a uniform external field can excite in a parabolic 
well is the well-known center-of-mass mo dJ, which involves no internal compression of the 
electron gas. We therefore focus in this paper on the electronic response in sing ldHand 
doubletSt^ square wells where no such restriction exists. 

The need for going beyond the ALDA becomes evident if one is interested in the linewidth 
of these electronic excitations. The ALDA accounts for Landau damping of collective modes, 
i.e. decay into single particle-hole pairs@. For high-frequency, long-wavelength modes, 
however, this decay mechanism is not effective, and damping is instead induced by dynamical 
xc effects such as multipair production. In other words, outside the regime where Landau 
damping occurs, the modes calculated within the ALDA will come out undamped. This 
limitation is overcome in non-adiabatic theories. In particular, in the formalism of Vignale, 
Ullrich and Conti (VUC)0, dynamical exchange and correlation lead to the appearance of 
viscoelastic stresses in the electron fluid, with complex and frequency-dependent viscosity 
coefficients depending on properties of the homogeneous electron gas. The viscosity causes 
an additional damping not contained in the ALDA. 

In an earlier paperc3, the VUC formalism was used in a model system of two-dimensional 
quantum strips, where it predicts a substantial inhomogeneity-induced enhancement of the 
damping of collective modes compared to plasmons in the two-dimensional electron gas. The 
interest in this model system is based on two reasons: first of all, it is analytically solvable, 
and second, it satisfies the conditions of validity of the approach, namely |Vno(r)|/n (r) <C 
kp(r) and uj/vf(t), where n (r), kp(r) and vp{r) are the local equilibrium density, Fermi 
momentum and velocity. However, no experimental results are available to compare the 
theoretical predictions with. 

The purpose of the present paper is to compare the performance of various non-adiabatic 
theories in computing the linewidth of electronic intersubband transitions in single and 
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double square quantum wells - systems in which high quality experimental data existllZrEI. 
Our goal is to clarify the limits of validity of the different approaches, i.e. the ones by Gross 
and KohnB and by VUC (the latter is equivalent to the method by Dobson et a/.0 in the 
situation considered here). For instance, the formalism of VUC is expected to be exact 
in the limit of slowly varying density (see above), and in particular does satisfy the HPT. 
However, this does not necessarily mean that it is the best method under all conditions: it 
will turn out in this paper that the VUC method works well only for those situations where 
the electronic motion is sufficiently collective. If, on the other hand, a single-particle-like 
behavior prevails (characterized by large gradients in the velocity field), then the method of 
Gross and KohnB gives results which appear to be in better agreement with experiments. 

The paper is organized as follows: In Sec. II we introduce the single and double quantum 
wells used in our comparative studies and calculate their electronic ground state. Sec. Ill 
contains a summary of the linear response formalism and the various approximations for 
the linearized xc potential. In Sec. IV we then calculate the lowest collective intersubband 
transitions for the single and double quantum wells and compare the results to experiment. 
Conclusions are given in Sec. V. 



II. ELECTRONIC GROUND STATE OF SINGLE AND DOUBLE SQUARE 

QUANTUM WELLS 

Before discussing the linear response in quantum wells, we need to find their electronic 
ground state. This problem has been dealt with by several authors beforeffi§'@S, and 
we give here a summary of the standard procedure to keep this paper self-contained. 

We work in the constant effective-mass approximation, taking a value of m* = 0.07m 
for the effective electron mass in GaAlAs. We furthermore employ a value of 13.0 for the 
dielectric constant e, so that the electronic charge becomes e* = e/^/e. In the following, we 
choose units such that e* = m* = H — 1. This means that lengths are measured in units of 
1 Oq = 98.3 A, and energies in units of 1 Hartree* = 11.27 meV. The difference between the 
lower conduction band edges of GaAs and Alo.3Gao.7As will be taken as 250 meV. 

We assume that the quantum wells have been grown along the z-direction and are trans- 
lationally invariant in the x — y plane. The Kohn-Sham eigenf unctions for an area A are 
then discrete in the z direction and plane wave-like (with wavevector q||) perpendicular to 
it. They can be written as 

^n^) = ^j^ Vl (z) , (1) 
with the Kohn-Sham eigenvalue 

E(q lbJ ) = ^+e J . (2) 

The orbitals (fj(z) satisfy the following Kohn-Sham equation: 

fid 2 \ 
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Here, the external potential v ext (z) is prescribed by the design of the quantum well. For the 
xc potential v xc (z) we use the LDA of Vosko, Wilk and NusairH The Hartree potential is 
obtained by direct integration of Poisson's equation asll 

z z' 

v H ( z ) = -4tt J dz' J dz" n {z") + 2nN s z . (4) 
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The sheet density N s is assumed to be a given constant (typically around 10 11 cm 2 ). To 
determine the density profile in ^-direction, no(z), we write 

n (z) = 2Y / ^(z)9(e F -E(q ll ,j)) 

= ^ E vp'Kti-^), (5) 

where e F and q F are the Fermi energy and momentum. The sheet density may thus be 
expressed as 



J dzn Q (z) = ^q F N< 




N s = dz n {z) = — q 2 F N occ - 2 V eA , (6) 



where N occ indicates the number of occupied orbitals, and the orbitals themselves are as- 
sumed to be normalized to one. To close the self-consistency cycle, we now have to determine 
N occ . To do this, we solve Equation (H) for q F : 
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1 v OCC 1 ' OCC j 

We now start with iV occ = 1 and, if necessary, keep increasing N occ until the condition 

is satisfied. In the cases we are interested in (see below), only the lowest subband is occupied, 
i.e. iV occ = 1. 

Let us now discuss specific examples of a single and a double square Quantum well, both 
of which have been experimentally studied. The single GaAs square wellS'El is taken to have 
a width of 384A. Note that this value is still within the ±4% range around the design-width 
400A of the sample used in the experiments^!. The well is filled with electrons up to a sheet 
density of N s = 0.97 x 10 11 cm -2 . In the following, we consider only the case of zero gate 
voltage and neglect any built-in static electric fields. 

The (asymmetric) double quantum well&B consists of two GaAs wells of widths 85A 
and 73A, separated by a 23A barrier of Alo.3Gao.7As. Its electronic sheet density is N s = 
2 x 10 11 cm -2 . Again, we consider only the case of zero electric field across the sample. 

We have solved the Kohn-Sham equation (|3|) for the two systems within the LDA. For 
our choices of N s , in both cases only the lowest level is occupied in the ground state. The 
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distribution of the bound levels is illustrated in Figure 1, which also indicates the total 
effective Kohn-Sham potential of the systems. The energy scale is chosen such that the 
Fermi level lies at zero. We see that the single square well supports nine bound levels whose 
spacing grows with increasing energy. The lowest subband spacing is obtained as E i2 = 8.18 
meV. The double well, on the other hand, is found to have only four bound levels. Here, 
the subband spacings are E\2 = 11.7 meV, E\z = 109.2 meV and En = 154.5 meV. 

In the following, we take the ground state of these systems and search for the lowest 
collective intersubband transition using linear response theory. 

III. LINEAR RESPONSE FORMALISM FOR QUANTUM WELLS 

A. Response equation 

The linear density response within TDDFT is given by 

n l( r, W ) = /^ K8 (r,r>)Mr>), (9) 

where XKs( r > r '> w ) is the noninteracting density-density response function, which in our case 
reads 

^k N , M (r)^ ; (r)^ |h/1 (r')^ kj| ,,(r') 



X KS (r,r» = 2 £ (/„ - /„) "' "' " . , (10) 

V w - [E(k\\,fi) - E{k\\,v)\ +irj 

k| | ,k(. 



with the usual Fermi occupation factors / M and /„, and 

v cS (r,u) = v eXk)1 (r,uj) + v s ,i{r,u) + f XCi i(r,^) . (11) 

Here, f ex t,i(r, w) is the frequency-dependent external perturbation, and v h< i(t,u) and 
f XCj i(r, uj) are the linearized Hartree and xc potential. Since our quantum wells are transla- 
tionally invariant in the x — y plane, we define 

ni(q||,z,w) = Jd 2 r\\ e^^n^Y.u) , (12) 

so that the response equation @ may be transformed into 

rii(q||,z,a;) = J dz' x KS (q||, z, z', u)v eS (q\\, z', a) . (13) 

In the following, we shall only be concerned with the case that q\\ = in Eq. (13). The 
linearized Hartree potential is then 

z z< 

v H<1 (z,u) = -4vr J dz' J dz"m(z",uj) . (14) 



-oo — oo 



The first-order xc potential will be discussed below. 

For the Kohn-Sham response function, we follow Refs. and obtain 
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-/V occ OO 

XKa(q\\, z,z',uj) = J2 ^ F i*>(q\\>u)<P l i(z)<PiM(*/)<p u ( z ) ( P''W) > ( 15 ) 

/1=1 u=l 



with the definition 



(2tt) 2 7 [q||k|| +a / ^(<?||) +o? + W7 q||k,| + a^(g||) - a; - J 

and 

<V(g||) = y - (e M - e v ) . (17) 
Setting gn = in Eqs. (|T6|) and (|TT| ) and performing the integration over fcn, we end up with 

F^ = 0, U ) = -^^\ ^— — + -) . (18) 

7r [ e„ - £ M + u + vq e v — £ M - u. >-ir}) 

This, together with Eq. (fL5|) , defines the response function x KS (qii = 0, z, z r , u). In the 
following, we will evaluate Xks summing only over the bound states <p v , i.e. the ones shown 
in Figure 1. This approximation is justified by the fact that the collective excitations we are 
studying are well below the continuum threshold for the quantum wells under consideration 
and thus involve only the lowest-lying states. 



B. Approximations for the xc potential 

In our case, the general expression for v xc> i depends only on the z coordinate, and is 
given by 

v xcA (z,u) = J dz' f xc (z,z',uj) ni(z',u;) , (19) 

where the xc kernel f xc (z,z',u) is a functional of the ground-state density n (z). 
The simplest possible approximation for t> xc i is the ALDA, which assumes 



J2 hom( \ 

f^ A (z,z',cu) = 6(z-z') - l " } 



dn 2 

n=no(z) 

where e^° m (n) is the xc energy density of the homogeneous electron gas, so that 

d 2 et m (n) 



(20) 



dn 2 



(21) 



n=n (z) 



In other words, the xc kernel in the ALDA has no frequency dependence at all and is purely 
real. This approximation is justified for cases where the external potential is slowly varying 
in time as well as in space. 

In a first step towards overcoming the restriction to slow variation in time, Gross and 
Kohni proposed the following plausible approximation (which we denote by DLDA, the D 
standing for "dynamic"), valid for the case where both n and n\ are slowly varying in space: 
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J xc,l 



(22) 



where f£° m is a property of the homogeneous electron gas, related to the dynamical lo- 
cal field factor as f^° m (n ,uj) = — lim q ^ 47cG(q,uj)/q 2 . Gross et al% have given a simple 
parametrization of f^° m (n , uj) (see also Ref. 0). A more elaborate calculation of f^° m (n , uj) 
has recently been performed by Nifosi, Conti and Tos@0 using an approximate decoupling 
scheme of the equation of motion for the current density, which accounts for processes of 
excitation of two electron-hole pairs. 

As already discussed in the introduction, the DLDA ( ^2f ) has recently been found to 
violate the HPT (for more details see Refs. jSHTop. However, a local frequency-dependent 
approximation which satisfies the HPT can be derived if the current density is taken as basic 
variable rather than the particle density. The resulting first-order xc vector potential in the 
form stated by VUC0 is given by 



^<Ur^) = VXcT(r^) 



r, u . 



(23) 



The dynamical correction to the ALDA xc potential is the divergence of the viscoelastic 
stress tensor 



a 



xc,i^u '/xc 



/ 2 



+ CxcV-u^. (24) 

Here, u(r,u) = ji(r, u)/no(r) is the velocity field, and rj xc (u , n (r)) and ( xc (u , n (r)) are 
complex viscosity coefficients whose explicit form is given in Ref. [H| 

In the one-dimensional case, it is straightforward to derive the associated xc potential©. 
It can be written as 



^f\ftrM I dz" n x {z» » 
n (z) J 



n Q {z ) 




f: 



dz" rii(z" , uj) . 



Here, we have defined 

tin 



2„homi 



d 2 e 



n 



dn 2 



n=n (z) 



(25) 



(26) 



is the DLDA ((22]) . The other terms are 



The first term on the right-hand side of Eq. ([ 
needed to satisfy both the HPT as well as Onsager's reciprocity theorem, i.e. symmetry 
under interchange of z and z' of the associated xc kernel fxc{ z , z',u) (see Ref. |]). 

Both the DLDA and the VUC approach lead to a complex xc potential. This means 
that the frequency of collective modes will come out complex as well, with a negative imag- 
inary part accounting for the damping. As a consequence, care has to be taken to perform 
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the proper analytic continuation of the Kohn-Sham response function and of f XC} both of 
which are analytic functions in the upper half complex plane, into the lower half complex 
planeEEi. 

The analytic continuation of the response function (|T5|) with is straightforward. 
This is because there are no branch cuts in the frequency region we consider, and the 
depolarization shift will move the mode frequency sufficiently far away from the poles in 
Ffj,v(q\\ — 0, uj). Thus, Xks(<Z|| = 0,z,z',uj) is analytic across the real axis in the frequency 
region of interest. 

For the analytic continuation of / xc , we proceed as shown Ref. The idea is the 
following: let us assume we have a parametrization of the imaginary part of f xc (uj) on the 
real frequency axis. Denote this parametrization by A xc (u), where v is real. Then, the 
expression for / xc (^) in the upper complex uj plane is 

oo 

/S^Woo + i f ^Idv, Zuj>0. (27) 

71 J V — UJ 

— oo 

Here, is the large-frequency limit of f xc . The analytic continuation into the lower complex 
uj plane is given by 

/ x » = foo + ~ 1 du + 2iA xc {cu) , %uj < . (28) 

7T J V — UJ 

— oo 

Note that the third term on the right-hand side implies that A xc has to be evaluated at 
the complex frequency uj, < 0. Therefore, in moving down from the real uj axis into 
the lower half complex plane, one has to make sure along the way that no branch cuts or 
poles of A xc come across. Due to its simple analytic structure, this is always the case for 
the parametrization of Gross et a/.0 The expression of Nifosi, Conti and TosiS0, however, 
has to be handled with some caution for large values of \^suj\. 



C. A simplified approach for the linewidth of collective modes 



The linewidth of collective electronic modes can be obtained in a simplified manner0il. 
The idea is to first calculate the modes using the ALDA. Outside the regime where Landau 
damping occurs, the modes will come out undamped, and the mode frequency Q is purely 
real. The damping due to dynamical xc effects is then added on perturbatively. For the one- 
dimensional case, the resulting linewidth (twice the imaginary part of the mode frequency) 
is given by 



\vtjdzut(z,n)£ z o-x 



J dz no(z)\ui(z, Q)\ 2 

From Eq. (^3f), the derivative of the stress tensor can be written as 

d d 

TrO- xczz (z, CI) = n —v£i A {z, fi) - n ^-v XCt i(z, Q.) . 
oz oz ' oz 



(29) 



(30) 



We therefore obtain the linewidth within the DLDA as 



S 



-pDLDA 



dji 
dz 



v 



and within VUC as 



-iVUC 



VL J dz no(z)\ui(z, 



Jdznl(z) ^V/xTKO)! 



(31) 



n J dz no(z)\ui(z, 



(32) 



Pdlda ^gpgn^g on gradient of the current density ji, whereas T™ depends on the 
gradient of the velocity field U\. These quantities are related to the density response rii via 
the linearized continuity equation 



IV. RESULTS AND DISCUSSION 
A. Numerical results 

Let us now apply the linear response formalism outlined above to the quantum wells 
whose electronic ground states we have calculated in Section II. The goal is to study the 
collective intersubband transitions between the first two subbands. Our main interest lies 
in the damping of these modes. 

The standard procedure is to set the external potential f ex t,i in the response equation 
(|13|) equal to zero. The resulting integral equation for ni(z, u) can then be viewed as a 
complex eigenvalue problem. The solutions with real eigenvalue 1 are the collective modes 
of the system. The real part of the mode frequency, Q = ^ftu, indicates the position of the 
mode, and its lifetime is given by r = r -1 , where T = 2Qu. 

Alternatively, one may consider the photoabsorption cross section 

a(u) = $fa(w) , (34) 

where uj is real, and 

2 r 

a{uj) = — — J dz znx(z,u) (35) 

is the dipole polarizability associated with an external potential of the form vi(z,u) = 
Ezcos{ut). We expect to find for our quantum wells that a{uj) is peaked at the frequency 
of the collective intersubband transition Q, with a Lorentzian profile of FWHM T. 

We have performed both kinds of calculations for the single and double quantum well. 
In addition, we compare with the linewidths T p following from the simplified approach 
presented in Sec. Ill C. Results obtained with the various approaches and approximations 
are summarized in Tables | and |I|. The first columns give the mode frequency Q obtained 
from ALDA calculations for the collective modes. In both cases, Q is already quite close 
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to the experimental result shown in the last column (within 4.2% for the single and within 
4.9% for the double square well). However, as mentioned before, the ALDA does not include 
damping of the modes. This means that u{uo) would consist of a delta peak at Vt. 

The second and third columns in Tables | and |I| give results for Q and T obtained 
within the DLDA, using the parametrizations of Gross/Kohn and Nifosi/Conti/Tosi for 
f^° m , respectively. Results within the VUC approximation are given in columns four and 
five, again for the two different parametrizations for f£° m . The first and third lines show 
Q and T as obtained from a solution of the response equation for the lowest eigenmode. In 
the second and fourth line, we show the position and width of the associated peak in the 
photoabsorption cross section a{u). Finally, the last line gives the linewidths T p as obtained 
from the formulas presented in Sec. Ill C. 

Let us first consider results for the single quantum well, see Table |. We note that for 
each choice of the xc potential there is an excellent agreement between the results for Q and 
T computed with the three different methods. The results for T p are extremely close to the 
linewidths obtained using the full calculation for the collective mode or the photoabsorption 
cross section, which demonstrates the reliability of the linewidth formulas fl3~T|) and (|32h. 

We find that the Gross/Kohn parametrization for f£° m shifts the mode frequency Q 
from its value in ALDA closer towards the experimental value. The agreement is fairly good 
within the DLDA, but the improvement is only small if the VUC approximation is used. 
The Nifosi/Conti/Tosi parametrization, on the other hand, leaves the mode frequency Q 
practically unchanged, within DLDA as well as VUC. 

Compared to the experimental value of 0.53 meV, the linewidth T comes out about 25% 
too high within DLDA, with only small differences between the two parametrizations for 
fx° m (0.683 and 0.655 meV, respectively). By contrast, the VUC linewidth is roughly a 
factor of five smaller that the experimental linewidth. It thus seems as if the DLDA yields 
a better description of the damping of the modes in the square well. 

One has to emphasize, however, that the experimental linewidth contains not only damp- 
ing effects due to electron-electron interaction, but also caused by additional mechanisms 
such as finite temperature, scattering off impurities, and fluctuations of the width of the 
wellS. Therefore, our calculations should provide a strictly lower limit to the measured 
linewidth. Clearly, this is the case if the VUC approach is used, whereas, under this aspect, 
the overestimation of the electronic xc damping within the DLDA is substantial. Note that 
this effect would be even more drastic in a parabolic well: in that case, VUC (which satisfies 
the HPT) would correctly give zero damping, whereas DLDA would result in an unphysical 
finite linewidth. 

For the double quantum well, however, the situation is different. From Table [TT], we see 
that the DLDA results are at least of similar, sometimes even better quality than for the 
single well. Again, the mode frequency Q agrees better with experiment if the Gross/Kohn 
parametrization for f£° m is used. In that case, the linewidth is slightly smaller than in the 
experiment (1.00 meV compared to 1.17 meV). With the Nifosi/Conti/Tosi parametrization, 
it becomes even smaller by almost a factor of three (0.403 meV). 

The VUC method, on the other hand, clearly fails to provide an adequate description 
for the double quantum well. Within the Gross/Kohn parametrization, the mode frequency 
Q drastically overshoots the experimental value by about 40%. The opposite happens with 
the Nifosi/Conti/Tosi parametrization, which shifts the ALDA value for Q down by more 
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than 1 meV. In both cases, the linewidth is dramatically too high. We also observe some 
differences between the results obtained from calculating the complex mode frequency and 
the photoabsorption cross section. Moreover, the linewidth formulas ( |3TD and ([32]) seem to 
break down here. This obvious failure of the VUC method for the double well calls for some 
explanation. We shall therefore examine the nature of the electron dynamics in the two 
quantum wells more closely. 

B. Physical picture: collective vs. single-particle 

To get some feeling for the underlying physics of the dynamical processes involved, it is 
helpful to look directly at the local behavior of various quantities characterizing the response 
of the electron liquid. 

Fig. H shows the ground-state density n Q (z) for the single square well and the density 
profile of its lowest collective mode, rii(z), together with the current density ji(z) and the 
velocity field u\(z). Calculations have been done within the ALDA. The other approxima- 
tions for the xc potential give very similar results. All quantities plotted in Fig. [2] vary quite 
smoothly across the well. Note that for the case of a parabolic well, uq{z) and ji{z) would 
have an identical shape, so that the velocity field U\(z) would be uniform. Here, however, 
Ui(z) has a shape similar to (though slightly smoother than) that of n (z) and ji(z), which 
may be viewed as an illustration for the non-parabolicity of the square well. 

The same quantities have been plotted for the double quantum well in Fig. |3]. Here, 
the ground-state density no(z) has two peaks, the lower one at the narrower well and the 
higher one at the wider well. The dip in n (z) is at the position of the barrier between the 
two wells. The density response 711(2) does not deviate too much from the one found for 
the single well, see Fig. |2|, apart for a slight hump at the barrier. The strongest differences, 
however, show up in the current density and the velocity field. We see from Fig. [3] that j\ (z) 
is very smoothly varying across the system, apparently ignoring the presence of the barrier 
between the wells. The velocity field Ui(z), however, is strongly peaked at the position of 
the barrier. 

From the examples shown in Figs. ^| and ^, we see that the most sensitive indicator for 
the physical nature, or rather the degree of "collectiveness" , of the motion of the electron 
liquid is the velocity field u\{z). Two extreme cases have to be distinguished: (i) Collective 
motion, as realized in a perfect way in a parabolic well. The velocity field is then only 
slowly varying. This is the case on which the HPT is based, (ii) Single-particle-like motion, 
occurring for instance in a tunneling process through a classically forbidden region, like the 
barrier between the two wells. Such a process is characterized by a pronounced peak in the 
velocity field, or, more generally speaking, by the presence of large gradients in u\(z). We 
are thus led to consider more closely the gradients of the various quantities involved. 

Note that the domain of validity for the VUC method is restricted by the conditionsHl§ 
k <C uj/vf, kp and q <C ui/vp, kp, where k~ l and q~ l are the characteristic length scales for 
variation of the external potential and equilibrium density, respectively, and kp and Vp are 
the local Fermi momentum and velocity. For the problems considered here (calculation of 
the photoabsorption cross section), the variation of the external potential is negligible, i.e. 
the condition on k is fulfilled. In turn, a measure for q is given by the quantity \n' (z)\/no(z). 
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The top parts of Figs. [| and || compare the local values of kp and uj/vf with \n' Q \/no 
for the single and double well. For the single well, the condition \n' \/no < u/vp,kp is 
reasonably well satisfied for the central part of the density distribution, but not towards the 
borders of the potential well. Note that a similar behavior would be found for a parabolic 
well. For the double well in Fig. [5], by contrast, the same condition is much more strongly 
violated, due to the fact that there are large density gradients around the region of the 
potential barrier in the middle. 

From the definition of the velocity (see Eq. fl33|)), we can write 

(36) 



n' {z) j[{z) u[(z) 



no(z) ji(z) u\{z) ' 

The two quantities on the right-hand side are plotted in the lower parts of Figs. f| and [5|, 
again comparing with the local kp. We see that for the single quantum well, \u[\/ui lies 
consistently below For the double well, however, the situation is opposite: \u[\/ui 

has sharp structures around the position of the barrier and is much greater than kp, whereas 
is relatively smooth over the whole interior of the system. Note, finally, that for the 
case of a parabolic quantum well we would have n' /n = since here the velocity 

gradient vanishes. 

We thus see immediately from Eqs. fl3"T| ) and ( |5^D that r vuc ^> r DLDA for the double well is 



explained by the large values of \u[\/ui. This suggests an additional, more refined local crite- 
rion for applicability of the VUC approach, requiring that | Vui(r)|/wi(r) <C | Vji(r)|/ji(r). 
This condition is closey related to the hydrodynamical picture inherent in VUC. It may be 
viewed as requiring a high degree of collectivity of the motion of the electron liquid in the 
system. This also suggests a pragmatic way to deal with the question of which method to 
choose for approximating the dynamic xc potential for an arbitrary system: use VUC only 
in those regions of space where the electronic motion is predominantly collective, i.e. where 
| Vu 1 (r)|/w 1 (r) < | Vji(r)|/j!(r). In those regions where the velocity gradients are large (i.e. 
the motion is more single-particle like), use the much more robust DLDA. Note that this 
prescription automatically satisfies the HPT. The so defined VUC-DLDA hybrid scheme 
yields a linewidth which always lies below the linewidths of the two "pure" schemes, or at 
most equals the smaller one. 

We have applied the VUC-DLDA hybrid scheme to our quantum wells. From Fig. |] it 
can be seen that for the single well the VUC-DLDA scheme coincides with the VUC method 
since here the velocity field is always smoother than the current density. In contrast, Fig. |5| 
for the double well shows that within VUC-DLDA the central part is treated with the DLDA, 



whereas the VUC is used for the outer region. Numerical results are shown in Table pi . 
Compared to DLDA (see Table [TI|), the resulting mode frequency Q is somewhat smaller if 
the Gross-Kohn parametrization is used, and almost unchanged with the Nifosi/Conti/Tosi 
parametrization. The effect on the linewidth is more pronounced: for both parametrizations, 
it is reduced to roughly half of its DLDA value, and is now substantially smaller than the 
experimental value, as it should. 

Fig. ^ illustrates the transition between the "collective" and the "single-particle-like" 
regime. We calculate the linewidth and lowest collective mode frequency (using ALDA 
and the linewidth formulas of Sec. Ill C) for our asymmetric double quantum well with 
variable barrier height, from its maximum value used so far, down to zero (which then 
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defines a single square well of width 181 A). All the while, the electronic sheet density is 
kept constant. The calculations have been done using the Gross/Kohn parametrization (the 
other parametrization yields similar results). We see that for decreasing barrier height, 
the VUC linewidth shrinks much more rapidly than the DLDA linewidth. The crossover 
occurs around 0.3 times the maximum height. This is the region where the mode frequency 
becomes comparable to the barrier height. The VUC-DLDA results, on the other hand, 
lie consistently below the DLDA results, and merge with VUC for small barrier heights, 
indicating a high degree of collectivity. 

V. CONCLUSION 

In this work, we have carried out a comparative study of two approaches to describe 
the linear electronic response of quantum wells within a density-functional framework. The 
two methods are the dynamical extension of the ALDA by Gross and Kohn (DLDA) and 
the recent implementation of time-dependent current density functional theory describing 
dynamical exchange and correlation in the language of conventional hydrodynamics (VUC). 
The latter approach is rigorously justified in the limit of slowly varying density, in the 
sense that it satisfies the harmonic potential theorem and other exact Ward identities and 
symmetries, which in turn are violated by the DLDA. A third approach, recently proposed by 
Dobson, Bunner and Gross^, turns out to coincide with VUC in the case of one-dimensional 
inhomogeneity studied here. 

Our studies show that the VUC functional should be applied in general only with a 
certain caveat, since there are physical situations of interest where it fails to work. Its 
domain of applicability was found to be closely linked to the structure of the velocity field 
of the electronic motion. If the latter is sufficiently smooth - which indicates a high degree 
of collective motion - then VUC leads to sensible results. On the other hand, if there are 
large gradients of the velocity field (as observed around the nonclassical barrier region in the 
double well), then VUC is obviously inadequate. By contrast, the DLDA approach was found 
to be much less sensitive to the nature of the electronic motion and led to a reasonably good 
description of the double well, although it clearly overdamps the predominantly collective 
motion in the single well. 

We finally proposed a pragmatic answer to the question of which functional to choose 
in a general case. We introduced a hybrid scheme using either VUC or DLDA in their 
respective spatial region of applicability. This provides a robust method to describe damping 
of electronic modes of either collective or single-particle-like nature, as shown for the case 
of a double well with variable barrier height. 
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FIGURES 



FIG. 1. Full lines: effective Kohn-Sham potential along the growth direction for the electronic 
ground state of the single and double square wells. The electronic sheet densities are 0.97x 10 11 cm -2 
and 2 x 10 11 cm -2 , respectively. The dashed lines indicate the bound single-particle levels. The 
energy scale is such that the Fermi energy lies at zero. Only the lowest level is occupied in both 
cases. The collective transitions we are interested in take place between the two lowest subbands. 

FIG. 2. Top: ground-state density profile no(z) for the single square well, see Fig. |]. Middle: 
density profile n\{z) of the lowest collective mode, calculated within the ALDA. Bottom: associated 
current density j\(z) and velocity profile u\{z). The latter has been scaled with a factor of 0.05. 

FIG. 3. Same as Fig. ^, for the double quantum well. 

FIG. 4. Top: rate of variation of the ground-state density, \n' (z)\/no(z), of the single square 
well, versus local Jcf and ujjvp. Bottom: rate of variation of associated current density j\ and 
velocity profile U\. 

FIG. 5. Same as Fig. ||, for the double quantum well. 

FIG. 6. Ratio of linewidth T p and lowest collective mode frequency for the double well with 
barrier height varying from zero to one, where one corresponds to the full height. Calculations are 
done with ALDA and the linewidth formulas of Sec. Ill C, using the VUC, DLDA and VUC-DLDA 
hybrid scheme, as indicated. 
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TABLES 



TABLE I. Lowest collective intersubband transition of the single GaAs/AlGaAs square well, 
calculated with various approaches for the linearized xc potential. The superscript 1 (2) stands for 
the parametrization by Gross/Kohn (Nifosi/Conti/Tosi) for f£° m . Q is the frequency of the mode, 
L its width. The subscripts indicate the computational scheme: search for the eigenmode of the 
system as pole in the complex w-plane (m), calculation of the photoabsorption cross section (er), 
and using the simplified approach from Sec. Ill C (p). All numbers are given in meV. 



ALDA DLDA 1 DLDA 2 VUC 1 VUC 2 Exp. 





10.25 10.63 


10.23 


10.31 


10.24 


10.7 




10.63 


10.24 


10.31 


10.24 




r 

J- m 


0.683 


0.655 


0.128 


0.104 


0.53 


Lfj 


0.683 


0.639 


0.128 


0.104 




r P 


0.686 


0.636 


0.130 


0.104 






TABLE II. Same 


as Table D, for the double square well. 








ALDA DLDA 1 


DLDA 2 


VUC 1 


VUC 2 


Exp. 




13.85 14.24 


13.88 


20.64 


12.55 


14.6 




14.24 


13.89 


20.23 


12.89 




r 

- 1 - m 


1.00 


0.403 


8.55 


4.15 


1.17 


Lcr 


1.00 


0.401 


8.52 


3.52 




r P 


0.988 


0.406 


11.07 


6.50 




TABLE III. Same as Table |Tl], for the double square 


well, using the VUC-DLDA scheme. 




VUC-DLDA 1 




VUC-DLDA 2 




Exp. 




14.07 




13.87 




14.6 


r 

- 1 - m 


0.620 




0.194 




1.17 


r P 


0.605 




0.192 
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Fig. 1 
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Fig. 3 




Fig. 4 
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Fig. 5 




Fig. 6 
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